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> : ABSTRACT 
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' Energetic particles that undergo strong pitch-angle scattering and diffuse through 

a plasma containing strong compressible MHD turbulence undergo diffusion in mo- 
mentum space with diffusion coefficient Dp. If the rms turbulent velocity is of 
order the Alfven speed va, the contribution to Dp from slow-mode eddies is ^ 
^ ; {2p^vp,/9l) [ln(/vA/£>|| ) +2y- 3], where / is the outer scale of the turbulence, 0.577 

is Euler's constant, and D|| is the spatial diffusion coefficient of energetic particles, 
^ which is assumed to satisfy D|| ^ Zva- The energy spectrum of accelerated particles is 

O ■ derived for this value of Dp, taking into account Coulomb losses and particle escape 

^ ■ from the acceleration region with an energy-independent escape time. 

I ■ Slow modes in the D|| ^ ZvA-limit are an unlikely explanation for electron accel- 

^ ■ eration in solar flares to energies of 10-100 keV, because for solar-flare conditions the 

^ ■ predicted acceleration times are too long and the predicted energy spectra are too hard. 

The acceleration mechanism discussed in this paper could in principle explain the rel- 
i. atively hard spectra of gyrosynchrotron-emitting electrons in the 100-5000 keV range, 

■ but only if D\\ ^ Zva for such particles. 



1. Introduction 

The first part of this paper treats stochastic particle acceleration by slow modes in compress- 
ible magnetohydrodynamic (MHD) turbulence assuming efficient pitch-angle scattering. The sec- 
ond part presents the energy spectra of accelerated particles taking into account Coulomb losses 
and assuming that particles escape the acceleration region on an energy-independent time scale. 
The third part of the paper shows that slow-modes in the efficient-pitch-angle- scattering limit are 
not a viable explanation for electron acceleration in solar flares in the 10-100 keV energy range, 
although in principle they could explain the electron spectra at 100-5000 keV inferred from mi- 
crowave gyrosynchrotron radiation. 
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Throughout the paper it is assumed that energetic particles propagate through a magnetized 
plasma containing both "large-scale" turbulent compressible motions that vary over distances 
much larger than the energetic -particle gyroradius p_^, and "small-scale" waves that vary over dis- 
tances ~ Pg. The small-scale waves scatter particles, causing diffusion in physical space with mean 
free path Xmfp as well as diffusion in momentum space. The large-scale motions also cause spatial 
and momentum diffusion. Only those large-scale plasma motions on scales S> Xmfp are considered. 
Particle acceleration by velocities on scales » has been considered previously by a number 
of authors [e.g., Kulsrud & Ferrari (1971) Bykov & Toptygin (1982, 1990, 1993), Ptuskin (1988), 
Dolginov & Silantev 1990, Katz & Stehlik (1991), Chandran & Maron (2003)]. In the present pa- 
per, the large-scale velocities are taken to be slow-mode eddies in strong anisotropic compressible 
MHD turbulence. Such eddies are elongated along the local magnetic field, and involve velocities 
directed primarily along the local magnetic field (Lithwick & Goldreich 2001, Cho & Lazarian 
2002). Part of the motivation for considering slow modes is that under certain conditions they 
can contribute significantly to momentum diffusion. An additional motivation is for completeness; 
by comparison the effects of fast modes and Alfven modes are relatively well understood (e.g., 
Jokipii 1966, Kulsrud & Pearce 1969, Berezinskii et al 1990, Miller et al 1997, Schlickeiser & 
Miller 1998, Chandran 2000, Schlickeiser 2002, Yan & Lazarian 2002). 

If the small-scale waves are Alfven waves, the transport equation describing the energetic- 
particle distribution function / is (Skilling 1975) 

a 



dt 



dp- 



.V/ = V.(Z.„bb.V/) + (V.w)f| + |, 



where p is momentum, b is the magnetic-field unit vector. 



W = U-|- ( -(1 -{J 



VAb 



(1) 



(2) 



V+-FV_ 

is an effective "wave-frame" velocity at which the cosmic rays are advected, u is the large-scale 
plasma velocity, va is the Alfven speed, v+ and v_ are, respectively, the pitch-angle scattering 
rates associated with waves traveling parallel to and anti-parallel to the magnetic field, p is the 
pitch-angle cosine, (. . .) indicates an average over p. 



\2(v+ 



+ V- 



is the diffusion coefficient for particle motion along the magnetic field. 



1 



v+v_ 



(3) 



(4) 



2 yv++v 

is the momentum diffusion coefficient associated with the small-scale waves, and % is the rel- 
ativistic Lorentz factor of an energetic particle. Equations (2) and (4) show that either w 7^ u or 
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Dpp ^ (Schlickeiser 2002). Nevertheless, for simplicity the momentum diffusion associated with 
small-scale waves is ignored (D^^ — > 0), and at the same time, inconsistently, it is assumed that 
w = u. This approximation has been standard in studies of particle acceleration by velocities on 
scales > Xmfp [Bykov & Toptygin (1982, 1990, 1993), Ptuskin (1988), Dolginov & Silantev 1990, 
Katz & Stehlik (1991)]; it reduces equation (1) to 

|^+u-V/ = V-(Z)||bb-V/) + (V.u)|^. (5) 

Equation (5) is used to calculate the momentum diffusion coefficient arising from the large- 
scale turbulence. It seems plausible that a nonzero Dpp would approximately add to the momentum 
diffusion coefficient associated with the large-scale turbulence calculated from equation (5), but 
future work is needed to verify this. The spatial diffusion caused by the large-scale turbulence is 
neglected. It is assumed that the rms large-scale turbulent velocity is comparable to va and that the 
ratio P of thermal to magnetic pressure is < 1. The power spectrum of the slow modes is assumed 
to follow the theory of Lithwick & Goldreich (2001), which has received support from direct 
numerical simulations (Cho & Lazarian 2002). It is also assumed that VAfimin ^ ^|| ^ vaK where 
^^min is the length along the magnetic field of the smallest slow-mode eddies under consideration 
(Jmin is either the length of the eddies at the dissipation scale or several times X^fp, whichever 
is larger), and / is the length of the largest (outer-scale) eddies. Since the time for an eddy of 
length A,|| (measured along the magnetic field) to be randomized is A,||/va (Lithwick & Goldreich 
2001), the inequality VA^^min <S ^|| <^ vaI implies that the largest eddies are randomized before 
a particle diffuses through them, and that a particle diffuses through eddies of length dmin before 
such eddies are randomized. The assumption D|| » VA^^min requires that the particle speed v satisfy 
V » VA since D|| ~ vA-mfp and dmm > Knfp- 

The procedure used to calculate the momentum diffusion coefficient Dp arising from the 
large-scale velocities is to average over the slow-mode fluctuations, neglect spatial variations in 
the averaged distribution function /o, and approximate the eddies on scales comparable to / as 
a uniform background field Bq. The eddies on scales <^ I have velocities <^ va and magnetic 
perturbations <^ Bq and can be treated using quasilinear theory. It is found that 



P^Dp^ ] , (6) 



with 

91 



ln|^).2,-3 



for VA^^min < ^|| < VaL (7) 



where Y~ 0.577 is Euler's constant. The corresponding acceleration time scale, p^/Dp, is of order 
the large-eddy turnover time //va, and depends only weakly through the ln(/vA/£)||) term on p 



-4- 



and particle species. An approximate version of equation (7) is obtained using phenomenological 
arguments in section 2.3. 

Section 3 presents the time-dependent spectrum of accelerated particles assuming an ad-hoc 
model of particle escape from the acceleration region with an energy-independent escape time. 
Section 4 gives the steady-state spectrum for electron acceleration taking into account Coulomb 
losses, again assuming an energy-independent escape time. Section 5 compares the predicted 
spectra and acceleration times to observations of electrons in solar flares. The results of the paper 
are summarized in section 6. Much of the notation used in the paper is defined in table 1. 



Notation 



Meaning 



Xmfp scattering mean free path of energetic particles 

^thermal Coulomb mean free path of thermal protons 

A,|| length of a turbulent eddy measured along the magnetic field 

A,_L width of a turbulent eddy measured across the magnetic field 

Vg energetic-particle gyroradius 

/ energetic -particle distribution function 

N{E) number of accelerated particles per unit energy 

B magnetic field 

length along magnetic field of eddies at the dissipation scale 

Jmin ^11 or several times Xmfp, whichever is larger 

/ stirring scale (outer scale) of turbulence 

VA Alfven speed 

Cs sound speed 

P ratio of thermal to magnetic pressure 

u turbulent velocity 

V energetic-particle velocity 

D|| diffusion coefficient for motion along magnetic field 

Dp momentum diffusion coefficient 

Tesc escape time scale 

Tacc acceleration time scale 

% '^acc/'^esc 

Y Euler's constant, ~ 0.577 



Table 1: Definitions. 
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2. Stochastic particle acceleration by slow-mode eddies in strong compressible MHD 

turbulence 

This section reviews properties of small- amplitude slow magnetosonic waves and compress- 
ible MHD turbulence and presents phenomenological and analytic derivations of the contribution 
to Dp from slow modes in compressible MHD turbulence. 



2.1. Properties of small-amplitude slow magnetosonic waves 

It is instructive to consider the properties of a small- amplitude slow magnetosonic wave prop- 
agating in a uniform background magnetic field Bq = Bqz with frequency co, wave vector k, and 
-C k. Such a wave varies much more rapidly across the magnetic field than along the magnetic 
field, like slow modes in compressible MHD turbulence (see section 2.2). The wave frequency, 
polarization, and velocity divergence when <g A: are given by (Lithwick & Goldreich 2001) 



to = 



CsVaI^z 



1 + 0(| 



(8) 



and 



k u = 



0, 



1 + 



klUii 



(9) 
(10) 

(11) 



where u is the fluctuating velocity, k is taken to lie in the xz-plane, and Cs is the sound speed. Equa- 
tions (9) and (10) imply that the slow- wave velocity is approximately aligned with the background 
magnetic field. The total pressure perturbation (magnetic plus thermal) vanishes to second order in 
k^/k (Lithwick & Goldreich 2001). An important property of slow waves with kz<^k when P < 1 
(and thus Cs < va) follows from equations (9) and (11): k ■ u ~ kz\u\. 
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2.2. Slow-mode eddies in strong anisotropic MHD turbulence 

The anisotropy of magnetohydrodynamic (MHD) turbulence has been studied by many au- 
thors.^ In this paper, it is assumed that the ambient plasma is stirred at a scale /, and that the rms 
turbulent velocity is comparable to va- In this case, the rms magnetic fluctuation at scale /, denoted 
Bi, is comparable to any mean magnetic field in the system. The stirring excites an inertial range of 
turbulent fluctuations extending from the large scale / to a much smaller dissipation scale. Within 
any box of dimension <g /, the fluctuations can be decomposed into the three MHD wave polariza- 
tions with respect to the average magnetic field direction within the box, Biocal- For example, when 
P <^ 1, the velocity fluctuations aligned with Biocai are associated with slow waves. The slow-wave 
fluctuations within the box can be thought of as a collection of nested eddies, where an eddy is 
simply a volume of some specified width measured across the magnetic field and length X\\ 
measured along the magnetic field, for which the velocity variation across the width of the eddy is 
comparable to the velocity variation along the length of the eddy. For values of X.^ in the inertial 
range, slow-mode eddies are elongated along Biocai> with 

X|| ~x5_/^/V3, (12) 
and the rms velocity variation across a slow-mode eddy is 

~ VA I — 1 (13) 

(Lithwick & Goldreich 2001, Cho & Lazarian 2002). Equation (12) means that the velocity varies 
more rapidly across Biocai than along Biocai- Suppose the vector separation between two points is r, 
with r in the inertial range of the turbulence, and let the rms velocity difference between the two 
points be 5u. If r ± Biocah equation (13) implies that 5m ~ v\{r/iy/^\ if r || Biocah equations (12) 
and (13) imply that 5m ~ va(?"/0^'^^- Slow-mode eddies are randomized in a time ~ ?1||/va due to 
the mixing of slow modes by Alfven-mode eddies (Lithwich & Goldreich 2001). 

The dissipation scale in the directions perpendicular to Biocah denoted d^, and the correspond- 
ing parallel scale = cff^l^l^, are given by (Lithwick & Goldreich 2001) 

di. ~ rg^themial for P 1, and (14) 

^^11 ~ ^thermal forp>l, (15) 

'See, e.g., Montgomery & Turner 1981, Shebalin et al. 1983, Higdon 1984, Higdon 1986, Oughton et al. 1994, 
Sridhar & Goldreich 1994, Goldreich & Sridhar 1995, Montgomery & Matthaeus 1995, Ghosh & Goldstein 1997, 
Goldreich & Sridhar 1997, Matthaeus et al. 1998, Spangler 1999, Bhattacharjee & Ng 2000, Cho & Vishniac 2000, 
Maron & Goldreich 2001, Milano et al 2001, Lithwick & Goldreich 2001, Cho & Lazarian 2002. 
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where Athermai is the coUisional mean free path of thermal ions, and ion-neutral friction and radia- 
tive cooling are ignored. Equation (14) is surprising in that linear slow magnetosonic waves are 
strongly damped in low-|3 plasmas on scales smaller than the coUisional mean free path: i.e., when 
(3^1, the damping time of slow modes with < ?ithermai is ^ ^n/^s. comparable to the linear 
wave period. The reason that turbulent slow-mode eddies can persist on parallel scales < ^thermal 
when P ~ Cg/v^ <^ 1 is that the cascade of slow mode eddies is controlled by the Alfven modes, 
and the cascade time X||/va is much less than the linear slow- wave damping time (Lithwick & 
Goldreich 2001). 



2.3. Phenomenological estimate of Dp 

Since only that part of the turbulence with X\\ ^ X^fp is considered, and since the momentum 
diffusion associated with small-scale waves is ignored, the time derivative of a particle's momen- 
tum induced by the large-scale velocity u is given by (Ptuskin 1988) 

dp V-M 

For P < 1 , the rms velocity divergence of slow-mode eddies of width 'k± satisfies (Lithwick & 
Goldreich 2001) 

(|V-u|t)^/^-^. (17) 
The rms contribution to dp/dt from eddies of length X|| is thus 

I) (18) 
3. All/ 



If a particle interacts coherently with an eddy of length X\\ for a time At, it incurs a random 
momentum increment of rms magnitude Ap ~ {dp/dt) At. The contribution to Dp from eddies 
of size is ~ {Ap)^/At. When JminVA ^ D|| <^ Iva, particles are confined within eddies with 
> D||/va throughout the time X||/va required for the eddies to be randomized in the turbulent 
flow. For such eddies. At ~ A,||/va. Each eddy size between D||/va and / makes a contribution 
to Dp of ~ P^va/ (9/). The contribution to Dp from all such eddies is 



^ 91 \ Dn 



forJn^n<(Z)||/vA)</, (19) 
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in approximate agreement with equation (30) below. The contribution to Dp from eddies with 
A, 1 1 < D||/vA can be neglected. 



2.4. Derivation of Dp in the quasilinear approximation. 

Particle acceleration by slow-mode eddies on scales <^ I is now treated analytically for (3 < 1. 
A scale /' is introduced with Amfp <^l' <^l, and an approximation is made in which eddies with /' < 
< / are replaced with a uniform magnetic field 5/z, so that Biocai is everywhere approximately 
along z. In terms of the Fourier transform of the velocity u(k,f), and in accord with equations (12) 
and (13), the spectrum of the (homogeneous) turbulence is taken to be (Goldreich & Sridhar 1995, 
Lithwick & Goldreich 2001, Cho & Lazarian 2002) 

(u(ki,?i)-u(k2,?2))=nki)^-^*''i-'2l5(ki+k2), (20) 

with 

for (/')-! <k±< dl^ [P(k) = otherwise], 

8{x)^e~\'\, (22) 

and 



p(k) = ^r'/3,-'o/3jy;^!) (21) 



y, = kTr^lW~'^, (23) 
where k± = ■^k^ + k'^. Given equation (17), it is assumed that 

(ki •u(ki,?i) k2-u(k2,r2)) 
= -kj^ P(ki)e-T^«l^i-^2l5(kj +k2). (24) 



^If slow-mode eddies were oscillatory, as opposed to repeatedly randomized, then eddies larger than 0|| /va would 
contribute less to Dp than eddies of length D|| /va, and the logarithmic term would disappear from equation (19). How- 
ever, although V • u is oscillatory within any infinitesimal fluid element (since the density does not change secularly in 
time), the average value of V • u over the volume of a slow-mode eddy is randomized in the time /va since the eddy 
is completely shredded and mixed with its neighbors in a time ~ ^||/va by Alfven modes. In addition, an energetic 
particle moves from one fluid element within an eddy of length A,|| to a new fluid element that is chosen randomly from 
the surrounding eddy- volume of width X± and length in a time / va due to cross-field diffusion and the turbulent 
mixing of fluid elements by Alfven modes. 
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For Biocai °= z, equation (5) becomes 



The distribution function is written as the sum of two parts, 

/ = /o + /i, 



(25) 



(26) 



where fo = (f), fi = f — (/), and angled brackets denote an average over an ensemble of real- 
izations of the turbulence. For simplicity, it is assumed that (m) = and V/o = 0. The ensemble 
average of equation (25) is then 



^ = -(u.V/i) + ^(V.u)-3^ 



(27) 



Subtracting equation (27) from equation (25) yields 



dt II dz^ 3 ' 



(28) 



in which products of fluctuating quantities have been neglected. Such products are small compared 
to \dfi/dt\ since the slow-mode velocity is almost along z, u <^ va for eddies on scales <^ I, 
and \dfi/dt\ ~ VA\dfi/dz\. Upon solving for /i and substituting into equation (27), one recovers 
equation (6) with 



(29) 



WhenJi„i„< (D||/va)</[|, 



Dr 



2pV 
91 



In 



+ 2y-3 



(30) 



where y ~ 0.577 is Euler's constant. An estimate of the contribution to Dp from all eddies, includ- 
ing those with /[, < A,|| < /, is obtained by taking /[, I in equation (30). 



3. Time-dependent energy spectrum of accelerated particles assuming an 
energy-independent escape time and no Coulomb losses 

In this section, the energy spectrum of accelerated particles is presented for p > po, where po 
corresponds to a particle velocity of a few times va- Coulomb collisions are neglected, and particle 
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escape from the acceleration region is modeled in an ad hoc fashion by adding a loss term to the 
right-hand side of equation (6), 



a/ 1 a 



dt dp 



a/ 



dp 



f 



(31) 



where the escape time Tesc is assumed to be momentum-independent. In equation (31) and through- 
out this section, the subscript on /o will be dropped, with the understanding that / now refers to 
the distribution function obtained after averaging over the turbulent fluctuations. Equation (31) is 
solved for p> po and t >0 subject to the boundary conditions 



fip = P0,t>0) = A 
f{p >po^t = 0) = 0, 
f{p^^,t>0) ^ 0. 



with A constant, 
and 



(32) 
(33) 
(34) 



Equation (32) would apply, at least approximately, if Coulomb collisions kept the distribution at 
low energies fixed and only a small minority of the particles were accelerated. It is assumed that 

(35) 



Dp = 



''&CC 



with Tacc independent of p. From equation (30) it can be seen that Tacc for acceleration by slow- 
mode eddies in compressible MHD turbulence is of order 1/va and independent of p and par- 
ticle species, aside from a possible weak momentum and species dependence arising from the 
ln(/||VA/£)||) term. The solution to equations (31) through (34) is (Kardashev 1962, Schlickeiser 
2002)3 



f 



A 
2 



+- 



P_ 
Po 

-) 

PoJ 



-(3+v/9+43c)/2 



1+erf 



(-3+V9+4x)/2 



■ln(p/po)+xV9 + 4x 
2v^ 

\ _ ( Hp/po)+xV9 + 4x 
V 2v^ 



where 



"^acc 
t^esc 



and 



is the error function. 



JO 



dy 



(36) 

(37) 
(38) 

(39) 



^See also Parker (1957), Ramaty (1979), SchUckeiser (1984), Ball et al (1992), and Park & Petrosian (1995). 
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Fig. 1. — The distribution function f ai t = Xacc/10 (solid line), t = Xacc/2 (dotted line), and 
t = 2Tacc (dashed line) for Xacc = 4Xesc> when Coulomb losses are neglected. 

The solution in equation (36) is illustrated in figure 1 for % = 4. The solid line corresponds 
to X = 0. 1, the dotted line to x = 0.5, and the dashed line to x = 2. Roughly speaking, / approaches 
the power law A{p/ po)~^^'^^'^'^^'^^/'^ for p — pQ up to p — pmax ~ poe'^^'^'^'^^, with / dropping 
off sharply for p > pmax- The width of the cutoff in ln(/7//7o)-space increases like For the 
case X = 4 plotted in figure 1, / ~ p""^ for p < pmax = Poe^'^- When Xacc > T;esc> Pmax ~ e^f/Vwc^^^ 
and the de-facto time scale for particle acceleration is ~ ^XescXacc- This time scale is short for 
small Xesc because the only particles to reach high energies when Xesc <S Xacc are a small fraction 
that happen to be accelerated much more rapidly than average. 

The energy spectrum N{E) is the number of accelerated particles per unit energy. In the 
non-relativistic limit, it is given by 



(40) 

p=y/2mE 



N(E) = 4%mpf{p) 

where m is the particle mass. In the non-relativistic limit, N{E) approaches the power law E""^ 
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for E < p^Jlm, with ^ = (1 + •v^9 + 4%)/4; equivalently, a specified power law N{E) oc 
corresponds to % = 4q^ — 2q — 2. In the ultra-relativistic limit, N(E) approaches the power law 
for E < pmaxC, with q = (— 1 + y/9 + 4%) / 4; equivalently, for a specified q, % — 4q^ + 2q — 2. If 
escape is neglected, then % = 0; this corresponds to ^ = 1 in the non-relativistic limit and q = 1/2 
in the ultra-relativistic limit. 



4. Steady-state energy spectrum of accelerated electrons with Coulomb losses and an 

energy-independent escape time 

In the presence of Coulomb losses, steady electron injection at energy Eq, and an energy- 
independent escape time Xesc, the energy spectrum of superthermal electrons N{E) satisfies the 
Fokker-Planck equation (Park & Petrosian 1995) 

dN d N 

= (DeN) - ^ {AeN) - — +Sd{E-Eo)®{t), (41) 

where in the non-relativistic limit 

D, = !^, (42) 

A,= '!^±ipD,)-?^, (43) 
dp p 

n is the electron density, S is a constant, and @{t) is the Heavyside function. Assuming 

Dp=^, (45) 

with Xacc independent of p, the steady-state (t oo) solution to equation (41) for £■ > £'o is [Park 
& Petrosian 1995, equation (71)]^ 

N{E) = cix^^/^M{a,b, |P|x-3/2), (46) 

where ci is a constant, 

E 

x= 2' (47) 



^Equation (46) of this paper follows from the results of Park & Petrosian (1995) upon noting two minor typos in 
that paper. In particular, the second line of equation (67) of that paper should read 5± = [(a + l)/2] ±jU, while the 
second half of equation (68) of that paper should read Xq = [{a + 1)/2]^ -|- 6 (Petrosian 2003). The correct versions of 
the equations were given by Park (1995). 
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and 



S=l+^, (49) 

P = ^. (51, 

Equation (46) is illustrated in figure 2 for the parameters Eq — 0.002 mgC^, n = 2 x 109 cm 
X — 10, and Xacc = 5 s. For x ^ P^^^, Coulomb losses become unimportant, and the spectrum 
obtains the same steady-state power law as in section 3. For x < |3^/^, Coulomb losses cause the 
spectrum to steepen relative to its high-energy power law (Hamilton & Petrosian 1992). 




Fig. 2. — The steady-state energy spectrum of accelerated electrons taking into account Coulomb 
losses and an energy-independent escape time, assuming Eq = 0.002 mgC^, n = 2x 10^ cm~^, 
X = 10, and Tacc = 5 s. 
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5. Application to electron acceleration in solar flares 

In this section it is shown that the acceleration mechanism discussed in this paper is probably 
not responsible for electron acceleration in solar flares in the 10-100 keV range. On the other 
hand, the predicted energy spectrum is roughly consistent with the energy spectra of electrons in 
the 100-5000 keV range, although it is not clear that <^ /va for such particles. 

A simple and fairly standard flare model is considered, as depicted in figure 3, with features 
drawn from the works of, e.g., Sturrock (1966), LaRosa et al (1994), Tsuneta (1996), Miller et al 
(1996), Tsuneta et al (1997), Miller et al (1997), and Aschwanden (2002). The flare arises from re- 
connection between open magnetic field lines above a pair of footpoints. The reconnection layer is 
assumed to be Petschek-like, so that the reconnection rate is sufficiently fast. Most of the magnetic 
energy is released outside (downstream) of the reconnection layer, where the magnetic tension 
associated with newly reconnected field lines accelerates plasma away from the reconnection site 
in a "sling-shot" action. It is assumed that turbulence is generated in this magnetized outflow 
with rms velocity ~ va, and that the turbulence is similar to strong homogeneous MHD turbulence 
(section 2.2). Electrons accelerated by the turbulence stream towards the solar surface, evaporating 
chromospheric plasma and emitting hard bremsstrahlung x-rays at the pair of footpoints that anchor 
the magnetic flux tube occupied by the energetic electrons. The evaporated chromospheric plasma 
rises, filling a magnetic-flux-tube loop with hot relatively dense plasma that emits soft x-rays. As 
additional field lines reconnect, the hard x-ray footpoints move progressively outward and the top 
of the soft-x-ray loop rises. Time-of-flight measurements typically place the electron acceleration 
region well above the tops of the soft x-ray loops and also above the hard x-ray emission that has 
been observed above the soft x-ray loops in a number of flares (Aschwanden 2002). 

The electron energy spectrum varies between flares and between the footpoints and above- 
the-looptop regions. If N{E) is taken to be a power-law E^^ in the 10-100 keV range, then a 
typical value of q is 4, although values between 1.5 and ~ 10 have been observed (Alexander 
& Metcalf 1997, Miller et al 1997, Petrosian 2002). Above ~ 100 keV, the electron spectrum 
flattens, with values of q in the range 2.0 — 2.7 for the three flares modeled by Kundu et al (2001). 
Time-of-flight measurements indicate that the electron acceleration time is very short. For the flare 
of December 15, 1991, 1932 UT, Aschwanden (2002) found that the time to energize electrons 
to 511 keV is <^ 0.36 s; Aschwanden (2002) pointed out that even the 0.18 s acceleration time 
to 511 keV found by Miller et al (1996) is too long to be consistent with the energy-dependent 
arrival times of electrons at the hard-x-ray footpoints in this flare. 

The start frequencies of bi-directional and type III radio bursts indicate that the electron den- 
sity in the acceleration region is in the range Ug = 0.6 — 10 x 10^ cm~^ (Aschwanden & Benz 
1997). Models of the gyrosynchrotron emission suggest that the magnetic field strength at the 
tops of flaring loops is ~ 300 G (Kundu et al 2001). The magnetic field in the acceleration region 
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\^^^^ J ^ chromosphere 

hard x-ray footpoints 
Fig. 3. — Reconnection model for a solar flare with a single pair of footpoints. 

is probably somewhat smaller. For typical parameters, it is assumed that = 2 x 10^ cm^-^ and 
B = 200 G in the acceleration region, yielding va ~ 10^ km/s. It is assumed that the outer scale of 
the turbulence is / = 5000 km, and that the region of greatest turbulence intensity is of size ~ /. 

It has been hypothesized that microturbulence in the acceleration region, e.g., whistler tur- 
bulence, causes Xmfp for energetic electrons to be much shorter than the Coulomb mean free 
path (Miller et al 1997). If D|| <^ /va as assumed in this paper, then the time for particles to 
diffuse out of the region of strongest turbulence, is much longer than the acceleration 

time, Tacc ~ l/vA- On the other hand, the plasma undergoes a bulk flow at speed ~ va away 
from the reconnection site that advects electrons out of the region of strongest turbulence on a 
time scale //va- For particles with D|| <^ /va, transport out of the turbulent region is dominated 
by this advection, and the escape time Tesc is then ~ //va ~ 0.5 s, independent of electron en- 
ergy (Blackman 2003). If JminVA ^|| ^ ^va, then particle acceleration by slow-modes is ap- 
proximately described by equation (30). Taking /^ ^ / in equation (30) and Z)|| = /va/10, one 
finds that Tacc = P^/Dp ~ 10//VA ~ 5 S, and % = Tacc/^esc = 10. 
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Coulomb losses are the dominant energy losses for electrons in solar flares with £ < 100 keV 
(Hamilton & Petrosian 1992). If escape is modeled with an ad-hoc loss term as in equation (41), 
then the energy spectrum of electrons accelerated by slow modes approaches the form depicted 
in figure 2. For energies above a few keV, Coulomb losses can be neglected and N{E) °^ E^^, 
with q = 2. Such a spectrum is significantly harder than typical solar-flare spectra in the 10- 
100 keV range, as discussed above. For energies above a few keV, the time dependent spectrum 
behaves approximately as in section 3, with N{E) oc £—2 fQj- ^ ^ Ej^ax, with E^^x — Pmax/^"^ °^ 
g(2t/Tacc)\/9+43c = gi4t/T:acc xhe time required for Emax to increase from 5 keV to 500 keV is 
then ~ O.SXacc ~ 1-5 s. This time is more than a factor of 10 larger than the energization time 
implied by time-of-flight measurements for the flare of December 15, 1991, 1932 UT (Aschwan- 
den 2002). The inability of the acceleration mechanism described in this paper to explain the 
observed steep electron spectra and short acceleration times makes it an unlikely explanation for 
electron acceleration in solar flares in the 10-100 keV range. 

The predicted spectra are roughly consistent with the energy spectra of electrons in the 100- 
5000 keV range inferred from microwave gyrosynchrotron emission (Kundu et al 2001). The 
electrons responsible for this emission may have considerably longer acceleration times than the 
electrons involved in x-ray bursts. The acceleration mechanism discussed in this paper could 
explain these high-energy electrons, but only if Z)|| -C Iva for electrons in this energy range. 

6. Summary 

The momentum diffusion coefficient Dp arising from slow modes in strong compressible 
anisotropic MHD turbulence is calculated for the case that the rms turbulent velocity is ~ va, 
pitch-angle scattering of energetic particles is very efficient (d mm VA <C i)|| <C /va), the energetic 
particle speed v is » va, and |3 < 1. It is found that 

forVA^jfmin<^|| <Va/, (52) 

where Y~ 0.577 is Euler's constant. Aside from a possible weak momentum dependence asso- 
ciated with the ln(/vA/i)||) term, Dp oc p^^ implying that the acceleration time Xacc — /Dp is 
independent of p. In addition, Xacc is of order //va and approximately independent of particle 
species. 

Slow modes in the D|| ^ ZvA-limit are an unlikely explanation for electron acceleration in 
solar flares to energies of 10-100 keV, because for solar-flare conditions the predicted acceleration 
times are too long and the predicted energy spectra are too hard. The acceleration mechanism 
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discussed in this paper could in principle explain the relatively hard spectra of gyrosynchrotron- 
emitting electrons in the 100-5000 keV range, but only if D|| <^ Iva for such particles. 
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